Laser-induced splittings in the nuclear magnetic resonance spectra of the rare gases 
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Circularly polarized laser field causes a shift in the nuclear magnetic resonance (NMR) spectra 
of all substances. The shift is proportional to the intensity of the laser beam and yields oppositely 
signed values for left- and right-circularly polarized light, CPL -/+, respectively. Rapid switching — 
in the NMR time scale — between CPL-(- and CPL- gives rise to a splitting of the NMR resonance 
lines. We present uncorrelated and correlated quadratic response calculations of the splitting per 
unit of beam intensity in the NMR spectra of ^^Ne, *^Kr, and ^^^Xe. We study both the regions far 
away from and near to optical resonance and predict off-resonance shifts of the order 0.01, 0.1, and 
1 X 10~® Hz for ^'^Ne, ^'^Kr, and ^'^^Xe, respectively, for a beam intensity of 10 W/cm^. Enhancement 
by several orders of magnitude is predicted as the beam frequency approaches resonance. Only then 
can the effect on guest ^'^^Xe atoms be potentially useful as a probe of the properties of the host 
material. 
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Nuclear magnetic resonance (NMR) spectroscopy has 
become one of the most successful techniques for the anal- 
ysis of molecular structure [l| . Irradiation by circularly 
polarized light (CPL) from a laser has been proposed as 
a potentially useful technique of enhancing the resolution 
of NMR spectra [21 • Early experiments searching for the 
effect, in a solution of chiral molecules, reported shifts of 
the order of 1 Hz in a 270 MHz NMR spectrum 3]. Much 
of that magnitude can, however, be explained by heat- 
ing effects and only a residual of at most 0.1 Hz could 
be interpreted as originating in other mechanisms [j|. 
Theory for the laser field-induced NMR shifts was pre- 
sented by Buckingham and Parlett [a> 13) in terms of 
the inverse Faraday effect Q, i.e., induced magnetiza- 
tion caused by the CPL beam propagating through any 
sample. Since the laser-induced shifts become propor- 
tional to the square of the electric field of the laser or, 
equivalently, the intensity of the beam 0, B lal ) initial 
proposals involved using high-intensity pulsed lasers to 
obtain shifts of the order of GHz |3|. This is, however, 
ruled out in NMR experiments, even leaving aside the 
heating problem, due to the long lifetime of the nuclear 
spin states. Hence, intensities of the order of tens of 
W/cm^, obtainable from continuous wave lasers, must 
be considered in estimates of the order of magnitude of 
the induced shifts. Left- and right-circularly polarized 
light cause opposite shifts. If the laser field is switched 
between the two modes rapidly in the NMR time scale, 
the spectral lines are doubled, with splitting proportional 
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to the beam intensity. 

NMR of guest rare gas atoms can be used as a sen- 
sitive probe of the microstructure of the surrounding 
medium 9] . In this work, we present ah initio electronic 
structure calculations of the laser-induced splitting per 
unit of beam intensity in the NMR spectra of free atomic 
^^Ne, ^"^Kr, and ^^^Xe, both at frequencies away from and 
near to optical resonance. In the off-resonance region, at 
standard laser wave lengths and intensities, the splittings 
in ^^^Xe are below the limit of current experimental ca- 
pabilities, while a large enhancement is expected when 
approaching resonance. 

A CPL beam of frequency u induces a magnetic field at 
the position of the nucleus K ina molecule given by 0,|a| 
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Here, E represents the electric field of the right (+) or 
left (— ) circularly polarized beam, and 6^^^ can be calcu- 
lated as a third-order perturbation expression from [lOJ 
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where we have used atomic units and the notation for the 
quadratic response function [ll| {{A;B,C))uii,uj2j equiv- 
alent to time-dependent third-order perturbation theory 
with the time-dependent operators B and C, as well as 
their respective frequencies wi and uj2- In Eq. Q, a is 
the fine structure constant, r^ and r^ are the compo- 
nents of the electric dipole moment, rix is the distance 
between the electron i and the nucleus K, and £iK.a is 
the a component of the angular momentum of the elec- 
tron i around the nucleus K. While —oP'Ylif-iK/'rlK i^ 
the magnetic field at nucleus K due to the orbital motion 
of the electrons, Eq. Q corresponds to the modification 



of this quantity by the external, time-dependent electric 
field. 

The relevant quantity for NMR experiments in the 
isotropic gas or liquid phase is the isotropic rotational 
average 
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where ea/j^, stands for the Levi-Civita tensor and {x, y, z) 
is the molecule-fixed Cartesian frame. For spherically 
symmetric systems, Eq. ^ reduces to one of the compo- 
nents of fo^«^, e.g. h^y^. The induced field is stationary 
and oriented along the direction of propagation of the 
beam. It couples to the magnetic moment jkIk of nu- 
cleus K , to give a term in the NMR spin Hamiltonian (in 
frequency units) in the high-field approximation as 
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where jk is the magnetogyric ratio of K, Ik.z is the 
component of Ik along the external magnetic field, and 
Eq is the amplitude of the electric field associated with 
the laser beam. Eq. Q corresponds to equally large but 
oppositely directed frequency shifts, A and —A, for the 
two differently polarized beams. A depends on the in- 
tensity Iq of the beam, and becomes l^l 
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where c and eo are the speed of light and permittivity of 
vacuum, respectively. 

Calculations of A/Iq were performed with the Dal- 
TON program [li|, using the implementations of analytic 
quadratic response functions of Ref. 14 at the ab initio 
self-consistent field (SCF) and multiconfiguration SCF 
(MCSCF) levels, and Ref. [I^ using density-functional 
theory (DFT). We studied the basis set convergence of 
the laser-induced shifts, starting with the uncontracted 
Gaussian basis set denoted by HIVu in Ref. 16]. Dif- 
fuse (small-exponent) functions were added to each block 
of angular momentum until convergence of the results 
at the correlated MCSCF level [using the complete ac- 
tive space (CAS) ^3 wave function specified in Table P 
to within 0.1% was reached. This procedure has been 
shown to lead to results close to the basis set limit in 
the calculation of nuclear magnetic shielding constants 
of the rare gases [l^ . The final, converged basis sets are 
(15sllp7d5/), (19sl6pl5d5/),and(26s20pl6d3/)for 
Ne, Kr and Xe, respectively. All the results reported 
below correspond to these converged basis. 

Both CAS and restricted active space (RAS) [l^ -type 
MCSCF wave functions were used. The details of the 
chosen active spaces are given in Table Q] The three- 
parameter hybrid B3LYP functional |23| was used in the 
DFT calculations. The calculations were performed for 
a set of typical laser frequencies. We also calculated the 
lowest singlet excitation energies for each atom with a 



CASSCF wave function, to estimate the range of fre- 
quencies where optical resonances occur. These ranges 
are 0.63-1.05 a.u. for Ne, 0.41-0.52 a.u. for Kr, and 0.35- 
0.45 a.u. for Xe. A/Iq was evaluated at frequencies ap- 
proaching the resonance, too. 

The results for A/Io in ^^Ne, ^^Kr and ^^^Xe are pre- 
sented in Table UTI The range of the shifts, calculated at 
SCF level, due to a laser beam intensity of 10 W/cm^ and 
wavelengths between 13190 and 4880 A, is 1.7-4.7 nHz 
for 2iNe, 9.6-30 nHz for ssRr, and 80-290 nHz for ^^SXe. 
The corresponding values for ^^^Xe are of the opposite 
sign and smaller than those for -'^^^Xe, and can be ob- 
tained from the latter by multiplying by the quotient 
of the respective magnetogyric ratios 7131x0/7129x0 — 
-0.296. 

Inclusion of electron correlation at the CAS level in- 
creases A/Iq by about about 15% for Ne, but decreases 
it by ca. 20% for Kr and Xe. Comparison of the results 
obtained with the different RAS and CAS wave functions 
shows that the choice of the active space is more impor- 
tant than using the multireference (CAS) wave functions. 
With RAS wave functions, correlation increases the shifts 
by up to around 15% for Ne and 25% for Kr and Xe, com- 
pared to the SCF level. For all atoms and frequencies, the 
DFT/B3LYP results are significantly larger than those 
corresponding to the ab initio methods. We have also 
tested the LDA and BLYP functionals, which resulted in 
even larger shifts. 

The larger shifts obtained for Xe as compared to Ne 
and Kr can be attributed to the larger polarizability of 
xenon. The electron cloud is deformed in response to 
the external electric field, thus giving a corresponding 
larger induced magnetic field at the nuclear site. Results 
of correlated CAS calculations of A/Iq and the dynamic 
polarizability a(uj), for Ne, Kr, and Xe as a function 
of the laser frequency uj, are depicted in Fig. ^ The 
frequency ranges from 0.1 a.u. to near the threshold of 
optical resonance, where the properties diverge because 
of the poles in linear and quadratic response functions. 
The growth of A/Iq with uj, although qualitatively simi- 
lar to a{uj), is faster than in the latter property. There is 
an enhancement of A/Iq by several orders of magnitude 
when optical resonance is approached. 

Finally, we have also carried out a preliminary inves- 
tigation of the infiuence of relativistic effects on A, by 
including the mass- velocity _ff™^ = —^ct'^J^i^i' ^^'^ 
Darwin iJ^ar _ 5 a^Tr^^ Zi<- ^- (5(rix) Hamiltonians 
as additional perturbations. The relativistic corrections 
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were evaluated using cubic response functions |2l| at 
the CAS level. The corresponding corrections to the 
shifts, A™^ and A^'^'^, have opposite signs and partially 
cancel each other. As expected, the relativistic effects 



in Ne are completely negligible. For °^Kr, b^^ and 
^Dar represent, individually, corrections of the same or- 
der of magnitude as the uncorrected value oibx, although 
their partial cancellation finally leads to values approxi- 
mately 25% larger than the nonrelativistic value. Finally, 
6x0 and fexc' are larger than the nonrelativistic values, 
roughly by a factor of three. Their cancellation leads to 
a corrected value of about —50 nHz for all the frequen- 
cies studied, with /q = 10 W/cm^. No definitive con- 
clusions may be drawn based on these results, however, 
as the response functions retain their nonrelativistic pole 
structure in this approach. Furthermore, picture change 
effects J23 on the hyperfine operator are also presently 
neglected. A fully relativistic quadratic response calcu- 
lation would be more appropriate, and will be pursued 
in the future. 

In summary, we have calculated the shifts induced by 
circularly polarized laser beam, to the NMR spectra for 
atomic ^^Ne, ^'^Kr, and ^^^Xe using first principles elec- 
tronic structure methods. At typical beam intensities 



and laser frequencies, the shifts are much too small for ob- 
servation, with the shift for ^^^Xe predicted at the order 
of magnitude of 1 /iHz. Experimental techniques exploit- 
ing higher beam frequencies would benefit from the dra- 
matic enhancement expected at near-resonant frequen- 
cies. If realised experimentally, the effect may provide 
a new characteristic signature of molecular structure in 
NMR spectroscopy. 
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TABLE I: Active atomic orbital spaces of the MCSCF wave 
functions used. 



Atom Wave Function A'^e" Active Space 
Ne CAS 8 2s2p -^ 3s3p3d 

RAS-I 8 2s2p ^ 3s3p3d 

RAS-III 10 Is2s2p -> 3s3p3d4s4p4d4f5s5p 

Kr 



Xe 



CAS 


8 


4s4p -^ 4d 


RAS-I 


8 


4s4p -^ 4d 


RAS-II 


18 


3d4s4p -^ 4d4f5s5p5d 


RAS-III 


26 


3s3p3d4s4p -> 4d4/5s5p5d5/6s6p 


CAS 


8 


5s5p -^ 5d 


RAS-I 


8 


5s5p — > 5d 


RAS-II 


18 


4d5s5p -^ 4f5d6s6p6d 


RAS-III 


26 


4s4p4d5s5p -^ 4f5d5f6s6p6d7s7p 



" A^e is the number of correlated electrons. 



TABLE II; Calculated laser-induced NMR line shifts per unit 
of laser beam intensity, A/Jq [in 10^^ IIz/(Wcm" 
gas atoms Ne, Kr, and Xe. 



Nucleus g;" (a.u.) SCF CAS RAS-I RAS-II RAS-III B3LYP 



^'Ne 0.0345439 


0.17 


0.19 


0.19 


- 


0.19 


0.31 


0.0428227 0.21 


0.24 


0.23 


- 


0.24 


0.38 


0.0656249 


0.32 


0.37 


0.36 


- 


0.37 


0.60 


0.0773571 


0.38 


0.44 


0.43 


- 


0.44 


0.71 


0.0885585 


0.44 


0.51 


0.50 


- 


0.51 


0.83 


0.0932147 0.47 


0.53 


0.52 


- 


0.53 


0.87 


^^Kr 0.0345439 


0.96 


0.80 


0.80 


0.85 


1.49 


1.76 


0.0428227 


1.17 


1.01 


1.01 


1.06 


1.86 


2.24 


0.0656249 


1.92 


1.60 


1.60 


1.76 


2.98 


3.62 


0.0773571 


2.34 


1.97 


1.97 


2.13 


3.57 


4.42 


0.0885585 


2.82 


2.34 


2.34 


2.56 


4.26 


5.27 


0.0932147 


2.98 


2.50 


2.50 


2.77 


4.53 


5.70 


i29Xe 0.0345439 


7.99 


5.32 


5.32 


9.58 


10.65 


19.17 


0.0428227 


10.11 


6.92 


6.92 


12.25 


13.31 


24.50 


0.0656249 


17.04 


12.25 


12.25 


20.24 


22.90 


41.00 


0.0773571 


21.83 


15.44 


15.44 


25.56 


28.22 


51.12 


0.0885585 


26.62 


19.17 


19.70 


31.42 


34.61 


62.84 


0.0932147 


29.29 


20.77 


21.30 


34.08 


37.81 


68.16 



" The frequencies correspond to wavelengths A — 13190, 
10640, 6943, 5890, 5145, and 4888 A, in the respective order. 



FIGURE CAPTIONS 



FIG. 1: Results of correlated CAS calculations of (a) the 
laser-induced NMR shift per unit of beam intensity A/ Jo, and 
(b) frequency-dependent polarizability a{uj), as a function of 
the frequency to, for atomic Ne, Kr, and Xe. 
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